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It lias teen shown that the hodograph eguations of motion can he 
derived in a symmetrical form, hy the choice of the velocity and the mass 
velocity as independent varlahles. The equations obtained hy the use 
of the velocity potential ^ the stream function, or their transforms as 
the unknown function are of the same general form and therefore can he 
treated in the same manner. 

Particular sets of solutions have heen studied independently of the 
gas law adopted and some properties of the series obtained hy TtiewuH of 
these sets have heen discussed. Approximate gas laws for which the 
solutions of the hodograph equations can he easily found have heen 
briefly discussed. 

The equations have heen further transformed so as to have as 
independent variables the complex velocity and the complex mass velocity. 
Two new generalized potential ftmctions can then he introduced that 
satisfy very compact equations. From these functions, all the 
quantities concerned with the representation of the motion can he 
derived hy means of formulas independent of the gas law adopted. By 
means of the generalized potential functions some developments have heen 
performed with the approximate Chaplygin-Ton K&rm£nr-Tsien law. 

An approximate transonic method has also heen suggested. 


INTRODUCTION 


From a purely mathematical point of view, the ordinary hodograph 
equations for the stream function or for the velocity— potential function 
and the eq'uatlons relating them to the physical coordinates are 
sufficient for the study of two-d3 menslonal isentropic flows. However, 
from a more physical point of view they are not very elegant because of 
their lack of symmetry in contrast with the symmetry of the corresponding 
relations for the incompresslhle case. 

^At present at Guggenheim Jet Fropulslon Center, School of 
Engineering, Princeton University, Princeton, New Jersey. 
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NoVj tlie eq^uations that defiixe the relcxjity potential 0 and the 
stream function -i}f are 

0n = O 

0s = V 
= 0 


^ = = m(v) 

where the suhscripts denote the differential quotients with respect to 
the element of streamline da or the element of normal dn, obtained 
from ds hy a coimter-clockwise rotation of 90 °, w, p(w), and 

Pq = p(0) represent, respectively, the velocity, density, anfl stagnation 
density. Equations (1), which are symmetrical with respect to ds and dn 
if p = Pq, conserve their property of symmetry for variable p if the 

mass velocity m is considered in some way the counterpart of- w. If 
the hodograph equations for 0 and ’i'’ (or for other functions) can be 
expressed so as to make w and m (instead of the relation connecting 
them) appear explicitly, the eqmtions will then have a symmetrical 
form that can be Interesting not only from a formal point of view but 
also from the fact that it can give rise to many possible developanents, 
some of which are Illustrated in the present paper. In particul^, it 
is possible to choose as new Independent variables the complex velocity 
and mass velocity and to lntrod\ice a new generalized potential function 
satisfying a veiy contact equation from which 0 and ijr, their Legendre 
transforms X and <d, and the physical coordinates x and y can be 
deduced by simple differentiations. 
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HDD0G3EIAPH EQUATIOHS 


Tlie hodograpli eoLuations can. "be directly deduced as follows. If N 
and S (fig. l) are the normal and the snhnormal to the streamline 
and 0 is the direction of motion at a point P 


z = X + iy = e^®(S + iU) 


dz = ei®[^ - W d0 + i(djr + S d0)^ ' 
= e^®(ds + i dn) = ^ 


(2a) 


(2h) 


where the defining eq_nations (l) have "been used in the last step. It 
follows from eq[uations (2) that 


d0 = w(dS - d0) 

dt= m(dIT + S d0) 


(3) 


These are two relations between esact differentials and therefore can he 
written as 


dS = dw + Sg d0 = i dw + 00 + njde 

dPr = %clm + n0d0=|-%^dm + (i^g- sjd0 


W 
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vliere m anfl y are two related variatles so tliat thie Tneanl n g of 
partial differentiation with respect to m is 


B ^ ^ ^ 

^aa dm. ^ m* ^ 


Since dS and dH are exact differentials , it follows from 
eqnations (it) that 


and 



so that 




> 




( 5 ) 


But from equations (it) 


mlV = ^ = ^w 


and 


’wSv = 0w 


or 





MCA TN 2lt-32 


Hence, putting 
transformation for 




eq^ual to 
results in 



and performing tHe same 


the following equations: 



1 St 




> 


which yield the we2JL-4oiown Chaplygin equations for 0 an<^ t in 
symmetrical form. 

It is seen from equations that 


and 



/ 


% 



and, with the aid of equations (5)> 


ow 


-Wg = mSj^ + S = 


SGns) 

cttn. 


and 
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Tliese eq.uations can be satisfied by putting, respectively. 


S=X^ , ^ =Xq , W = c%i ^ mS = -<1:0 (7) 

Eq^uations (7) are consistent if 


= -wo^ > -0^ = ni)C^ 


(8) 


a symmetrical system of equations in X and oi very similar to 
equations (6).^ 

From equations (3) and ( 7 ) it is deduced that 
d0 = w dXy. — Xg d9 = d(wX^ — 'X) 


and 


dilf = m a<i^ — CC 0 dS = d(ma:^ — cd) 


Hence, to id-thin an unessential constant. 


0 = »X„ - X 

/ (9) 


■»lr = mo:^ - oj J 

3 

wbicb give 0 and ijr in terms of X and co. Tbe functions X 
and <D are of course tbe Legendre transforms of 0 and ■'1*' considered 
as functions of tbe physical coordinates. 


Tbe functions X and cu are distinguisbed by tbe fact that once a 
solution of equations (8) is Imown all tbe other functions concerning 

^Equations (8) have already been written in tbe present form by 
Bateman (reference l) . 

3Relatlons (9) already have been derived in tbe present form by 
Bateipan and Peres (references 1 and 2), 
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the phyalcal representation of niotion can he derived hy simple 
differentiations, 0 and ilr helng obtained from eq^uations (9) and z 
from equations (2) and (7). However, If 0 and ijr are the known 
functions satisfying eq^uatlons (6), integrations are necessary to deduce 
the other q.uantities. After Integration and determination of the 
constants so as to satisfy eq_uation (8), the following e3p)licit 
expressions for X and cd are obtained from eq[uationa ( 9 )j 


0(w^,0l)sin(9 - 0i)d0i 


+ C^ cos 0 + Cq sin 0 



£D = -m 



ilfCwi,0)d 





ijf(wj,,0-]^)sln(0 — 02^)d0i 




0(w^,0l)cos(0 


0jJd0^ — sin 0 + Cg cos 0 


where Wj, and 0j- are two arbitrary reference quantities, = mCwjv), 

and C-^ and C2 are two arbitrary constants with no influence on 0 

and ilf. It is readily deduced with the aid of equations ( 2 ) and (7) 
that 



- 1 


“r 


w„ 


e^^^d©^ - (^Ci + iCgj 


By differentiation. 


dz = e 


10 








which agrees with equation (2b). 
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Rrom equations (6) and (8) tlie hodograph. eq^uations for the 
functions Xj and ijr can. "be easily shown to he as follows: 



Each of these equations induces to the Laplace equation if m = w. 

The following eq.mtions are obtained from, equations (10 ) with 
as independent variahle: 




(i + m2) 

Sw2 


■ ^ 

+ (l 


: 0 


8w 





(i- 

• m2) 


- 

II 

o 

(i- 

■ M^j 

W M + (l , 

- 1^) 

^=0 

U 

8w 

80 ^ 


(l 

II 

%l 

1 

0 


8w 

80 ^ 




(10a) 

(lOh) 

(10c) 

(lOd) 

w 
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wliere M, the Mach nTnnberj is defined "by (see eq_"uation. (15)) 


1 _ m2 = 

d log V 


For every partlcxilar law m(w); that is, for every p(w), ezjlicit 
equations are obtained. For M <1 (subsonic flow) the equations 
the elliptic type; for M > 1 (supersonic flow) the eg.uations are of the 
hyperbolic type-. 

The hodograph equations are frequently transformed so as to 

dm ~ 

simplify the second-order terms. .Thus, if for ^ 


= (v - [” ^(s) “ ° '' ^ 
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and eo[uations (lO) are transformed Into 


+ he ~ ^ ^ ° 


(13a) 


+ %e * - 0 


(13t) 


■ d log a 
+ ^0 + ° 


(13c) 


d log a 

+ ^e ~d5C — " ° 


(13d) 


If 0 = P0^j ^ ” F ^ “ a and cb = acojf.^ the equations 
for 0^j ■'I'ljj., and oq^ are readily shown to he 




hxk + ^^*00 ~ ^ ° 




^*00 


. 49 — 


^XX a , 2 ” ° 


dX 


°^XX + “*00 “ ^ “2 

dX 


{l) = ° 

/ 


(14) 


For — < 0 analogous transformations can he performed with the 
dw 

introduction of dX-j^, and defined in the same way as dXj a, 

and 3 with -dm. instead of dm. It is seen that X represents some 
"kind of an integral mean between log w and log m. Introduc ing the 
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/^p^l/2 

sound velocity a = {-^) and the Mach number M = ^ 

\dfy ^ 

following eq[uations deduced from the Bernoulli eq^uation 


leads to the 
dp + pw dw = 0: 


X 

a2 + w2 = 0 

d log p 


^ ^og P ^ > 

d log w 


ii2S^ = l_M2 

d log W / 


(15) 


Then for M 1, 


, = \/l — M^ d log w 

(I6a) 


(l6b) 


( 16 c) 


There is obtained a corresponding set of eq[uationB for dX^j and 
when M > 1 hy simply replacing 1 — hy — 1 . 


Approximate Methods 

It is seen from equations (13) that for constant a or 3 these 
ecLuatlons reduce to the Laplace eq,uation. The first possibility is to 
be rejected for subsonic motion because, as equation (l6) shows, it 
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giTes p as an Increasing function of M. Tlie second one is th.e well- 
knovn Eirm^n-Tsien approximation (references 3 anil 4) 


= fO >4-1^ = ^^K 

P 


(17) 


"wiiere K is a constant. This eq^uation reduces to tlie Chaplygin 
approximation for K = 1 (reference 5 anil derived studies). 

Some considerations that will he useful in a subsequent section are 
now introduced. With the help of equations (12), equation (I 7 ) can he 
immediately integrated to obtain 



^ = k(i + HwS) 


where H is an arbitrary constant. Integrating now equation (l 6 a) anil 
the Berno u l li . eqTiations (l5)j there follows, respectively. 


X = log 


w 


1 + 


x: 


+ Constant = —log 


Hw 



+ Constant (19a) 


and 



2^ fO 

EH p 


Constant 


( 19 b) 
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I^om eq^uations (l5)j for a real isentropic gas. 



and. 


d. / ^ ^ _1_ 

\p^/ a^ 


( 201 ) 


where 



la the TTW.TrtTmnn 


a^ is the stagnation sound, velocity, and 


velocity at zero density, 

7 is the adiabatic index. 


A comparison between eq^uations (l8) and (20) is shown in figure 2 
where the law, equation (l8), is represented by a strai^t line, 
Chaplygin takes it as the tangent to the graph of equation (20) 

at w = 0, so that K = 1 and H = In. the K^tnn^r-Tsien method^ 

the tangent is taken at w = the speed at infinity, so that 


2 



y 


^This method has been often presented in a less coherent form, as 
the constants of equations (l8) and (I9) are determined for different 
conditions, thou^ the formula for the correction of pressure coefficients 
is not affected by this incoherence. 
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01)serTe that for K 1 the telItis of p at w = 0 differs from the 
exact stagnation density pQ. It Is easily seen that the value for KH 

In eq^uatlons (21) gives the slope of the true Isentroplc relation at a 
point corresponding to the conditions at Infinity, hut, as this depends 
only on the value of KH, It Is seen also that HiCririaa*s condition Is 
satisfied for every parallel to the tangent at v = (fig. 2), that 

Is, for the KH value given hy equation (2l) hut for different values 
of K. This suggests the possihility of improving the mean approximation 
of the Kjarm^nr-Tslen method hy an appropriate choice of K. The KarmSn- 
Tsien correction formula for the pressure coefficient has then to he 
modified. The loodifled formula. Independent of the constant of 
integration in eq^uation (l9a), is then 



•which reduces to the Karm^a-Tsien correction formula for the value of E 
given hy eq.uation (21). Application of the modified formula -with some 


value of K he-tween 


and 


1 


gives -values of 


Cp in 


better agreement -with experlmen-fcal values. 


For sui>ersonlc motion the hodograph equations reduce to the 
simplest hyperbolic equation (-wave eq-uation) for ap = Constant 

or 3q = Constant -with cip and Pp given hy equations (12) -with —dm 
instead of dm or hy equation (l6) ■with — 1 in place of 1 — M?. 
The second possihility is now to he rejected because the resulting -value 
of p increases -with M. The first possibility gives 
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that is, after integration.'' 


= E2. (% ~ "w^) 



v2 = ^^^0^ 


The constants and can he detennined so as to satisfy E&rm^’s 

condition: 


% ■ 2 =^ (>^.2 - 1 ) 


■'0 


%% = a^M 


PO 




mJ-1 


/ 


With = Constant, eq.mtions (l3c) and. (I3d) (modified for M >1) 

reduce to the simple vaTe eq[uation. The general solution can therefore 
he represented hy, for instance. 


X = f(Xi + 0) 


( 22 ) 


5 Peres has already indicated a law of this kind (reference 2 ). 
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Other laws for which the hodograph equations reduce to the Laplace 
equation for M < 1 (or to the wave equations for M > l) eire imme- 
diately deduced from equations (l4), (or from the corresponding equation 
for the supersonic case) as the laws which make one of the four quant- 
tities a, l/a, 3, and l/p, (or those corresponding for M > 1) linear 
in X (or in X^). For the true isentropic gas the curves of 3, 

1/ttj and 1/3 as functions of X are shown in figure 3 . For w = 0, 

X = — ooj and a = 3 = 1 the Chaplygin and the Karmhn-gslen approximations 
replace the true shape hy a horizontal line which can give an approx- 
imation not too had even for w = 0. Approxi imte laws for which 
equations (l4) reduce to the Laplace equation are represented hy 
arhitraiy (and generally not horizontal) strai^t lines. It is seen 
that for X = — o> these lines diverge hopelessly from the tme law. 

Hence these laws do not appear to he convenient for the approximate 
representation In a large range of velocity, nevertheless they can 
possibly have application when the variations of velocity from a mean 
value (for instance , the value at infinity) are small. In this case it 
is possible to achieve a better approximation than with the Karmknr- 
Tsien method hy taking as the approximate law the tangent to one of the 
curves of figure 3. Hie resulting approximate p^p curve will have a 
contact of second order with the real isentropic. 

For the supersonic case^ a]_, 3i, 1 /° 1 j functions 

of X 2 _, are shown in figure k for 7 = l.^l-. An interesting possibility 
is given by the ciirve of which can be well approximated by a 

strai^t line between M = 2 and M = 10. Hence in this range the 
exact hodograph. equation in ~ differs very little from the 

simple wave equation in X^^ and 9. Therefore the general splution of 

the supersonic motion in the said range of M is approximately 


X = ^ f(Xi + 9) 


where ct^ and X^^ are given by the true isentropic law. This 

approximation seems to be better than the approximation given by 
equation ( 22 ) . 
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Exact Solutions of The Hodograph. Eq^uations 

•povrer set .— Manj studies hare teen developed, on two sets of 
paxticular solutions of the todograpit eq^uation in 0 or ''I'’ ^ the so- 
called power set and exponential set^ characterized ty the fact that 
for the incompressitle case they reduce respectively to the natural 
powers of the logarithms of the complex velocity and those of the 
complex velocity itself. The symmetrical form of eq^uations (10) makes 
it possible to present these solutions in a form that seems interesting. 

By the introduction (thou^ it is not strictly necessary) of four 
auxiliary q_uantities ZSj X, and satisfying in the same order 
eq^uations (10), and some supplementary conditions, the four complex 
OLuantities 

F = X+ IS 


F = X + icD 


G = 0 + it 

and 

3 = 0 + it 


can he made to satisfy not only eq,uations (lO) in the same order tut 
aJLso the relatione 


and 


mF^ = iF0 


wF„ = iFfl 

m 0 


(23a) 

(23t) 


h - -“e 


( 23 d) 
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correspanding to eq^uations (6)^ (8), and the relations 


G- = •wl' . — F 

W 


G- = _ F 


J 


(2k) 


corresponding to eq_nations ( 9 )* Conversely F and F are given ty 
fonnnlas corresponding to those written in section entitled "Hodograph. 
Equations" for X a nd , cd: 


F = -w 



G(w^,0l)sin^0 - 0i)d0i 


+ 


“r 



G(wj._,0i)cos( 0 — 0]^)d0i 


( 25 ) 


pnd a similar eq_uation for F with. m. and w, G and G interchanged. 

The q.nantities x j ^ determined so that they reduce 

to X, (D, and i}r for the incon^aressihle case m = w; hence the 
eq.uations 


Fi = F^ = )^ + lo>i 


Gj^ = GjL - 0i + 


and 
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satisfy tlie Laplace eq^uatlons 


and 


= 0 


= 0 


( ■where A = ^ | to which ec[mtions (lO) reduce ■when m = w. 

V S(log w)^ he^J 

An operator (( )) is now defined as that "which, when applied to 

/ vh 

flog — I (w_ heing an arbitrary reference velocity), transforms it 
into 

(ih 0 - ^ /c ? /r ? • • • ‘ 



formulas ( 26 c) and ( 26 d) could preferably be "written as operations 

and log 


on powers of log 






and are different from the 


operations given by eqmtions ( 26 a) and ( 26 b) on powers of —log — 


w^ 


and —log 


m 

n^* 
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I'rom theBe equations there is deduced 







and the analogous equations obtained hy interchanging w arn^ m, and/or 
hy changing v into l/w and m into l/m. 


It is then easily verified that, with Iq = 1, each of the following 
functions 


F = (((i°8 ^ ^ 

= “)“i) 

O = (((l°8 ^ 


7 

is a solution of the corresponding eq^uation, ' the signs having been 
selected so as to satisfy also eq^uatlons ( 23 ). The choice of the upper 
or lower sign does not affect the values of anfl oo derived 

from eq^uations ( 28 ) so that for the solution of the flow problem it is 


7 

'This kind of solution has been first discovered by Bergman (see 
for instance reference 6 ) and by Bers and Gelbart (reference 7 ). The 
present form is new and more symmetrical. 
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sufficient to retain, for instance, only the upper onejS nerertheless in 
some cases it can he useful to consider solutions with hath signs. 

The solution represented hy the values, eq^uations (28), of G 
pnfl S does not coincide with the one corresponding to the values, 
equations (28), of F and for it differs from the values of G- 
pnfl 8 f derived from F an<^ S' hy means of equations (2k) . This can 
he shown as follows: From the recurrence formula, which is easy to 

verify. 



and from the values, directly deduced with Iq = 1 , 

- H ^ 

and 

l2(w,m) - I Ij^(m,w) = 12(14) + ^ Il(;4) 


There is deduced with the aid of equations ( 27 ) 

(1 - w - m Ih-l(^:»'«’) 



(29) 


and, the analogous relations obtained after interchanging the variables. 

75 " — — 

In this case for m = w the four soltitions, equations ( 28 ), reduce 

(to within a multiplicative constant) to (los ^ 10 ) , that is to the 


power set for the inconqKressible case. 
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Application of egnations (24) to eq^uations (2,8) yields tlie 
following eq^uation: 



and a similar eq.uation for -G-, if w is interchanged with, m, with, 
the ± sign. This G differs from the elementary solution, 
equations ( 28 ), althou^ it is a linear coaribination of such elementary 
solutions, a fact holding also for the inccmipresslhle case. In that 
case, however, the espression for G contains only an elementary 
solution with the upper siga if F is so, for then w^ = m^. 

The physical coordinates are easily deduced; for hy eq^uatlons ( 7 ) 


S + Ilf = Fj^ = 


i 

m 



> 


S + 


ixr = F^ = - F„ 
m V 6 


(30) 


where S and IT are auxiliary q^uantities connected with x 
hy relations similar to equations ( 7 ). Hence 


z = e^®(S + ilT) 


with 
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It is now again Terified that, as already ohserved, the choice of 
the sign does not affect the results concerning the values of 0, ijr, 
and z hut only the introduced auxiliary functions. In suhseguent 
work, therefore, only the upper sign is retained. 

Linear combinations, and in some cases infinite series of the 
elementary solutions, eguations (28), are still solutions of the 
corresponding equations. 


Infinite series in tho nover set .— If, say, F is given hy an 
infinite series, then developing and inverting the order of summations 
gives 



n 

^ Ili(w,m) 
h=0 


(n - h) ». 


h=0 n=i 


nl 

(n - h)l 


On(-i0)“^ 



Ij^(w,m) 


,d^7^=-i0 


( 31 ) 


where 


A(5) = 

0 


represents the function corresponding to this power series. 

Uow, independently of the convergence of this power series it is 
readily verified hy differentiation, that eguation (31), whenever it 
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convergsa, is a solution of the eq^uation in F (see equation (lOa)). 
Similarly when it converges , the series 


F = 


h=0 


Ih(m,w) 


[ d ^ V ^=-10 


(32) 


ohtained from equation (3l)^ty simply interchanging 
satisfies the equation in F (see equations (lOh)). 
there are the analogous solutions 





w and m „ 
For Gr and G 


(33) 


NoWj let G and G he deduced from F and F hy means of 
equations (2h)j with the help of eqmtion (29)* There results the 
following expression for G; 



and a similar expression for G after w and m, w^ and m^ have 
heen Intesrchanged. 

The values of 0 and ilf deduced from these expressions for G 
and are the same as those derived from equations (33) iT 


1 . 1 /^r “r\ dim 


B( D = -A(-n 
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The formula for I" can he -written In a different form: 
S’ = 2 i\(v,m) + ia^^^^Ce)! 

h=0 ^ J 

■witli 

A(-10) = ajL(e) + 18^(0) 

and 

B(i0) = -b^(e) + ihgCe) 

It follows that 
with 

hj^(0) = -a^(0) _ ^a 2 *( 0 ) 

and 

hgC©) = - 02 ( 0 ) + — a^*(0) 

Wp 

r\^ 

There are similar eq.tiationB for F and G. 
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Hence for w = 


\ 

= ai(e). 

> m 

(Or = 



J 

BO tliat A and »1 1 ^116 solutions of eq.uations ( 31 ) aiid- (32) are 
detennined "by the Talues of X and its radial derlTative on the 
circle w = w^. of the holograph plane ctr hy the corresponding values 

for CQ. Similar statements hold for B, (f>y and ijr and solutions of 
eq.uations (33) . Hence the solutions 'written depend upon "two arhi'trary 
functions and^ in their region of convergence, represent the general 
solution of the holograph eq,uations. Haturally they do not give any 
indication of the behaviour of the corresponding solutions at v = 0; 

for as m-> 0, “ as (log w(^ so that "the origin is 

certainly outside the region of convergence. 

In fact, if the solutions must he regular at v = 0, only one of 
the functions a(0) and h(0) can he chosen arhitrarily, and other 
representations of the solutions are needed to determine "the other. The 
region of convergence of the series, eq.uations (3l)j ( 32 ), and (33) ^ 
depends on -the form of A(^) and B(^). Hovever, a general idea of 
its shape can he given hy mak-lng very general assumptions about these 
functions and m(v) . 

Let r(0) he less than 'the distance in -the ^-plane between the 
point —10 and the nearest singulari'ty of A(^) and let A^y,y (0) 
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"be tlie upper iDoimcL of |a(^ )| on flie circle "witli center —10 and. 
radius r. Ttien tlie Cauchy’s Ineguality glYSs 


d^A(^ ) 




5 111 

^ = -10 


(35) 


■and a similar expression for B_, -with and r^^ in place of 

Ajjjg^ and r. For tiie if 

(k-J^ ] =(z] =f^2) 

^m. d log V/max '™'^iax ' ^ /u 


= a 


'max 


1 dm 


■w d log w 


= |m‘( 

max mpLT 




d(l/w) 

d ioftr/vT 


V 


d(l/m) 


d log (1/w) 


max 


■ (v)inax ° ' ° 


> (36) 


= d 


y 


are the upper hounds of the 'wri'tten gixantities between w and Wj., 
eguatlons ( 26 ) gire for even Talues of h 


|ll^(w,m) I < ^ |log 


Z_ 


|%(I'b)I < |i°e ^ 


( 37 ) 


and the same limitations for Iii(m,-w) and Iji(l/m,l/w) . 
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For odd values of h. tlie iq)per "bounds for 


and 






are obtained from the cccrresp ending espressions In 

equations ( 3 T) ty multiplying them, with (a/b)^^^ and (h/a)^^^, 
respectively^ and those for |lj^(l/w,l/m)j and |ljj^(l/m.jl/w)J are 

obtained hy multiplying the expressions in equations (3?) 
and (d/c)^^^j respectively . 

Consideration of the series (31) shows that the terms of the series 


k=0 




d2kA(? ) 


d^2k 


00 

^ = -10 k =0 


d2k+l^(^ ) 


cl£ 


2k+l 


^ = -10 


are less than the corresponding tenns of the geometric series 



which converges when the ratio is less than. 1. 

Hence the series ( 31 ) and (as it may be deduced in the same way) 
(32) converge absolutely in the region 


(ab)V2 jiog < r(0) (38a) 

Similarly, the series ( 33 ) converge absolutely for 

(cd)^/^ |log < rj^(0) (38b) 

How since ab and cd are functions of w and v^, and since r 
and are quantities idiich Increase with the distance fPcan the 

singularities of A and B, the general shape of the region of 
convergence in the holograph plane is a curved strip, which contains the 
circle w = w and whose width will be a minimum when —10 is nearest 
to a singularity of A(? ) or B(^ ) . 
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If a singularity lies on th© Imaginary axis, tliat is, if A(— i0) or 
B(i0), and their derivates have a slngular*ity for some ralue of 0, the 
corresponding width of the region of comrergence will he zero. This 
happens, for instance, when the reference Telocity is the Telocity 

at infinity of the flow round a hody. 


An ohserration of some interest is that as, for plausible laws, 
p/pQ is 0 ( 1 ) for M = 1 then (ah)^/^ and (cd)^^ (equations (36)) 


are o( (l - M^j 
and near 1. 


or 


0^|l 


(the larger of the two) for M 


Hence equations (38) show that for giTen A(^) and B (0 the 
width of the region of conrergence is the greatest and the rapidity of 
conTergence the best near the sonic line M = 1. Therefore it is 
belisTed that the solution represented by series (31) ^ (32), and (33) 
may haTe applications in the solution of transonic problems, naturally 
in combination with other methods conTerglng in the rest of the field 
of motion. 


Finally obserre that the deTelopment of equation 
handled differently so as to obtain a power series in 


( 31 ) can be 
0; 


^ n*cxjj^li3^(w,m) 

h=0 h! n?sh 


which by use of equation (27) and similar eapressions and with 
f^(w) = X n!o^Ij^(w,m) and f 2 (w) = 2. n!o^Ij^(m,w) becomes 


ru=0 


F = 


k=0 


n=0 

-sa 

> 

■k=0 


^ (x, M1\(w) + ^ V ^ fi ^fw MTf2(^) 

(21c)! [_ dm \ dw/J 1 ' ■'^q (2k +1)! dm dwV dm/J ^ 


Similarly, 


~ ^ ^ t ^ (v fe)Tfp{w) ^ m ^ fv ^(m |-Tl\(w) 

^ (2k)! L d-w \ dm/J 2^ (2k + 1)1 dw L dw/J 1^ 


The functions F and F, wheneTer the series conTerge, satisfy the 
corresponding equations and the relations ( 23 ), as can be directly 
Terified, with arbitrary fq and f 2 * Analogous solutions hold for G 

and G with 10 in place of — i0, l/m and l/w in place of w and m, 
fmd two arbitrary functions gj^(w) and g 2 (w). For real Talues 
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of f jl and. f£ tlie real and; imaginary parts of F and. are 
obtained, directiy; so that tlie eq.uation 


CO 




+ 


CO 


(-i 0 ) 


2 k+l 


kz =0 ( 2 k + l)I 



depends only on f 2 (w); cotrrespondingTyj and g^ 15000106 real and 
the equation 



depends only on gj^(w) . 

The physical coordinates are then found hy means of equaticoos (2) 
and ( 7 ) to he 


z 


2; + iy = e^® 


fy (-19)^ 
(a)i 



^ {2t + 1)'. I dv|_ dw^ j 
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From X + 1^ Idiere follows "by TtiRnns of eiuatiaa (9) an expression 
for 0 + lT}r wMch, must coincide with the written one if 


gl(w) = ^w 


dw 



or 




-w 


d(l/w) 


I^caa this coincidence the following interesting farnmlas can he deduced: 



and 


1 d 
V d(l/m) 




which are easy to yerlfy directly. 


The meaning of the written solution is readily found hy ohserv±ng 
that T^en the series converge, for 9 = 0 , then X = fj^(w), 0 = g^Cw), 


^1 • P P 

0 ) = ijf = 0, 2 = and y = 0, so that f^ =yx dw, g^ = Jv dx, and thei 

■vdiole solution is determined when the “axial” law of dlstrihutlon of 
velocity is gl-ven. 


Hence the solution, under a somewhat different and more explicit 
form, reduces to the one studied hy Li^thlll (reference 8) in hie 
work on the transonic flow in symmetrical channels . As Li^thill 
observed the coefficients of the expansions hecome infinite at sonic 
speed (for then dm = O), so that the series diverge in the transonic 
region. In this region, however, the solution can he found hy following 
the Lighthill's ingenious method, that is, inverting the series 
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giving lijr or iih .0 on© giying sine© ■fcli© co©fficl©irbs of "tli© 

(10) 

inT©rt©d s©ri©s are finit© at sonic spsod. Tt© application of Li^tMU’s 
metliocL can t© mad.© ©asi©r ty tli© pr©s©nt form of th.© solution. 

Til© ©xponontial s©t .- If in ©g.uations (3l), (32), and ( 33 ) 

A(£) = %a(“)(£) = 


and 

B(£) = 


■with, artitrary n^ tli©n tli© functions 




F = E^(w,m)© 
F = E^(m,-w)©‘ 


-inS 


-inj9 


ine 


ln 0 




y 


•with. co©fflci©nts d©fin©d hy 




^ 3 ^) = i n^Ij^(v,m) 


h =0 


Eji(m,w) = 5Z n^Iii(ia,w) 
h =0 


(39) 


m 
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33 


and. similar eq.uations for tlie other will "be solutions of the 

correspoad±Qg eq.uations (lO) and ( 23 ) "idilch reduce^ for m = w 

to (we~^®/wj.j (eq.uationB (lOa), (lOh)^ (53a), and (23h)) and 

to (equations (lOc), (lOd), (23c), and ( 23 d)), 

that is, to the esjonential set for the incan^essihle case. 

It is immediately verified hy means of equation ( 37 ) that the 
series converge for all values of w and for which a, h, c, 

nnd d are limited, that is, for which p and M are limited and not 
zero. For plausible gas laws this excludes only the values w = 0 

It is therefore seen, and easily verified directly, that 
eq.uatlons (4o) and the two other eguations for E^^ are, in order, 

solutions of the ordinary differential eq.uations 




1 

d 

|i ^ ' 

0 
II 

1 

(hlc) 

m d(l/w) 

(w d(l/la)j 

1 

d 

1 Hi 

0 

II 

4 

1 

(hid) 

•w d(l/m) 

[f d(l/w) J 


deduced from eq.izations (lO) hy taldng X, o>, 0, and as the product 

of a sinixsoidal factor in n9 and of the corresponding function 

or of w. For the normal isentropic law, equations (hi) 

become the known eq.uations of the bypergeometrlc type and have been the 
object of the investigations of many authors. (See, for instance, 
references 5, 6 , 8 , 9 , 10, 11, and 12.) Equation (hid) has generally been 
studied with particular regard to those solutions that satisfy the 

condition that is unity at w = 0 . 
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It Is iramsdiately verified that the solutions do not satisfy 

this condition and that the carrespanding canditions at v = are-* 



The symmetrical form, of equations (4l) allows some general 
relations to he easily derived. Some of these relationSj that may he 
useful for further develoigmenta, are now stated briefly. 


The can he considered as the supei^osition of two independent 

solutions^ of equations (4l): 


n 


(w,m) 

k=0 


k=0 


( 42 ) 


^Solutions of the kind of 
hy Bers and Gelbart, reference 


7 . 


and 0^ 


were first introduced 
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for wMch. the conditions at v = Vj. are: 


(CJr = 1 



(Sn)r = 0 


/ dSn(w,m) \ / dSn(m^ 

\ dw 4 ~ \ dm 


=n 

f-r* 


They are connected hy relations similar to those connecting the 
exponentials and the hyjerholic cosine and sine: 


C_^ = C 
-n n 


^-n “ ~®n 


± 8 

±n n n 


2C^ IP Eq + E 


-n 
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~ Sn(w,m)Sji(m,w) = 1 
I^Q.(w_,m)E_ii(m,w) + En(m,w)E_5i(w_,m) = 2 


m 


dCn(w,m) ^ , 

— to ~ ' 


m 


dSjj^(-w,m) 


dw 


= nCjj(mjw) 


m 


= nE^(m,v) 


Similar relations hold idien w and m are replaced hy -their Inverse 
values . 

All the solutions of eq.mtlans (4l) can he represented hy linear 
camhlnatlons of Cjj^ and hut these functions (as all the series 

In Iji) are not suited to give the hehavlour of solutions near w = 0. 
Since this behaviour Is very Important for many physical applications, 
it is necessaiy to follow a different method of Investigation: let 


and 


G = 
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"be solutions of the corresponding eguatlons (lO) -which reduce to -the 
exponential set for the Incompresslhle case. Then ^n^-’ 

and are solutions of equations (4l) . UoWj If for negative values 

of n the last two must coincide -with and ijr (such -that, for 

Instance, =1 at w = 0) and If for positive values of n the 

first two must coincide -with -the corresponding and cd^. It Is seen 

that Xq/w, Xj^/m, mT^^, and must he equal to unity at w = 0. It 

Is now shown that this Is possible for all values of n except some 
exceptional values : Let 


d log Zq 

®n d.w 


/V. d log i_ 

^ = ^ ' dm" 


^ w cL(l/m) 


? = 1 <3- log Xn 


m 


d(l/w) 


m 


It Is Immediately seen from equations ( 23 ) -that must satisfy 

Elccatl*s equation 


1 ^ 
n dw 


dm^ Tj 2 


( 44 ) 
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and that and Tjj^ must satisfy the cotrrespondliig equations 

■with m and w intercha n ged or inverted. Itirtheniiorej these 
qxiantitles are connected the following relations: 






(45a) 






(45h) 


nm — E_v 

T = I 

-n nwKji^ — m 


(45c) 


It is deduced from eq.uations (43) that 


w 



(46a) 



(46h) 


= exg 


= e2j) 



(46c) 


(46d) 


so that 
eOl the 


if 



these q.Toantities are eq.ual to unity and analytic at w = 0 if 
integrands are analytic there. It can he shown in fact that. 


P 


is an analytical function of 


w near 


V = 0, the integrands 


of eq.-uations (46a) and (46h) are zero nnfl analytic at v = 0 for all 
values of n except negative Integral and half— Integral values (for 
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nnTj r tile negatlYe integers if 


_P_ 

PO 


is an analytic function of 


•w^); 


wMle tlie integrands of eq.uations (l|^c) and (ll^d) are zero and analytic 
at V = 0 for all Talues of n except positive integral and ha.Tf— 
integral values greater than 1 (for only the positive integers greater 

than Ij if ^ is analytic in w^). 


These results reduce to the vell-ino'wn results when the eq.uations 
•are hypergeometrlc. In this case T.lght. hm (reference 8) has given the 
most ccm^jlete discussion of the solutions of the eq.uation in and 

deduced i inp ortant theorems, some of which may possihly he generalized 
following the present method. It can he seen that the exclusion of the 
pole at n = — 1 for the eq.uatlons in and is a general 

property, Tdiich does not hold for the eq.uatlons in and £%. 

For n = 1 the solutions of eq.uation (lf4) (and of the analogous 
eq.uatlons) and the corresponding solutions of equations (2k) are: 







F = = we"i® 




F = Z^e“^® = me"^® 


o ” = 5 




>19 = i .IS 


a = 1,6-^'' = — e- 
1 w 


(^Ta) 

(l^Th) 

(4?c) 

(4Td) 

(^Te) 

(kjf) 
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It is seen frcam. equations (24) that the values of G and G (hence 
of 0 and t}) corresponding to equations (47a) arwl (47h) are 
Identically zero. In fact, it is seen from, eqmtions (2) and (3) that 
this solution represents merely a displacement of the origin of the 
physical coordinates. Eq.uatlons (47e) and (47f) are more interesting 
as they coincide with the well-known Ringlet solution (reference I3) . 
The corresponding values of F and F are determined hy eq.uatlon (25) 
(and the analogous equations) . Thus^ 



plias a constant multiple of we and me The coefficients 

of e^® in these formulasj together with w and m respectively, 
represent two Independent solutions of equations (4la) nnt^ (4lb) 

for n^ = 1. The operators ^ ^ ‘which generally 

allow the deduction of two independent solutions of the equations 
in 0jj^ and from two independent solutions of the equations 

in Xn and oJq, suffer an exception for n = 1 as they produce only 
one set of solutions; that is, l/m and l/w, respectively; for when 
applied to w and m the result is zero. This exceptional case is 
explained in the section entitled "A new set." When n = 1 

and ^ equation (45c) ‘becomes indeterminate , hut the corresponding 

value of T__^ can he dedxiced from this relation as the limiting value 

for n 1 of the indeterminate expression. It follows from 
equation (44) that with the condition = 1 at w = 0 


En = nj 


^ miwi2(n l)dwj^ + o[(n - l)^ 
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4l 


Hencej It is deduced tiiat 


T-1 


Um 

n— >1 — m 



mjL dVjL 


wj da^ 


(49) 


and 




dnij 


Similar esjressions hold for T_^ a-nfl 

(r_i)-i = ii<i - I 

= +1 = I 

represent the second solution .of equations (4lc) and (4ld), independent 
of Yq = ^ and “ v and reducing to zero at w = 0. The 

espresslons 0je~^® and represent a kind of motion hetveen 

two paral l el walls. These solutions could he directly obtained hy 
inverting and exchanging the variables in equations (48) and 
putting Wj, =•%" = 0* 


T__ 2 . Hence the equations 


mn. 

Jo 


(50) 


£ 


™i 
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If n-^oo, eq.Tiation (44) bAows tliat 

d-w2 


Similarly 


tlie equations 


~ 2 _ 1 _ <i^ 

“ 2 gjij2 

J^OO 



(51) 


? 2 ^ _ 1 _= d(l/M.^) 

T_«? d(l/w2) 

J 

coincide -wltli l/a^ and l/P^ (equations (12)) and can te explicitly- 
calculated ly equation (l6). Then eq-uations (46) show -that 


2;*, -X. =Y 


1 = T 1 = e^- 


(52) 


•sdiere \ = log w — /;f‘ — (d log m/d log w)^/^ d log w coincides 

•with the value deduced hy integrating -the dX given hy equation (11) 

X 

and detennlning the constant of integration so -bhat — =1 at w = 0. 

V 

Hence e^ is -the suhsonlc asymptotic value of and 


for n-^ 00 . 
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The solirbions of eq.uatians (4l) jiist discussed are hourid to the 
solutions, equations (42), hj siiig)le relations. For Instance, 


Xn(~^) 


= Cn(w,m) + Eii(wp)Sii(w,m) 




ikn^'Wr) 




(53) 


J 


Hence as w 0 

and |Cjj^j and 

1 lTn 

Cn(w,m) 

_ .lini 

Sn(ni,w) 

0 

Sn(w,m) 

w-> 0 

Cji(m,-w) 

lim. 

"=( 14 ) 

lira. 


V-? 0 

^ 14 ) 

■W-5» 0 



= Vvr) 


The following Interesting expansions are deduced hy applying 
equations (53) to the solutions in closed form obtained for n = ±1 


V = VpC2(w,m) + m^Sj^(v,m) 


J 


(54) 


(and the corresponding eajansions "with the variables interchanged or 
inverted) . 
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Finally, it should, he ohserved. that differentiating (for instance) 
eq.uation (4la) ‘with respect to v yieldn 


1 ^ 
m d.w 


V 


d 

d(l/m) 



(n2 - 1) 


dX 


n 


dw 


= 0 


This equation is of the same general form as equations (4l) with 

2 2 

only one of the variahles inveirbed and 1 — n in place of n , it can 
therefore he treated in the same way as equations (4l) . 

Hence two particular solutions similar to eq.uatlons (42), 



can he defined throu^ the integrals Ij^ given hy the formula ( 26 ) 

hy sln^jly replacing the present variahles, Fccr n^ >1 the second 
series is imaginary and must he divided hy i to obtain a real 
solution. The general solution for d)^dw is given hy a linear 

combination of these solutions, and the general solution for is 

obtained by integrating and adding an appro xi mate constant. Now, this 
must coincide with the one in terms of Cjj^(w,m) and Sj^(w,m) . It is 

then easily derived that 


C 





+ 


•Wr 5n(^j~w) 
m n 
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and 



1 




n 


Analogous relations ottalned "by interchanging and inverting the 
variahles also hold. For n“^ = 1 these relations give eq.uations (5^) 
as a particular case. 

Other Interesting relatione can he deduced in the same way. 

Infinite series in the exponential set .— The series in 

and have heen used (as series in ■ and 

could he) hy many authors in the case of the normal isentroplc law 
(references 5j 9, and 12) . 

They seem to have their natural field of application in the 
problem of two-dimensional gas Jets, as Chaplygin first showed in his 
classical memoir. 

The apjlication to flows around bodies seems to he more 
difficult, especially for flows with circulation. The difficulty 
arises first from the presence In the holograph plane of a singularity 
at w = Woo from the ensuing necessily of enploying more than one 

series development in the exponential set with different seq.uences 
of n (as appears already in the incon5)reBsihle case) with added 
eventual terms in other sets, and of insuring that the different series 
axe the continuations of each other. 10 This can he achieved (although 
in a not Yevj single way) hy putting the condition of continuity of the 
solutions n-ntl of their derivatives on the the transition curves (often 
circles), as has heen done hy Tsien and Euo (reference 12) and as the 
author himself has done in an unpublished warh in a scmiewhat different 
way, hut it is believed that the main obstacle to this method arises 
from the difficulty of insuring that the body will have a closed 
contour idien a circulation is present. In fact if the so-called 
"natural" series (that is a series having the same coefficients as in 
a chosen incoirqpressihle case with circulation) is used in one part of 


lOrhis difficulty is avoided in the method by Bergman (reference 6), 
which uses a different type of expansion and uses the. series in the 
exponential set only as eventual auxiliary series. 
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tlie hodograpli plane (so tliat a "basic series from -wiilch the coefficients 
of the other series -will he deduced hy the foresald method is obtained), 
the resulting body -will be closed only in the limiting Incongpressible 
cane. In a tentatlTe method the author has tried to obtain the closing- 
up of the contour by talcing the coefficients of the basic series as 
stable fuuctions of (or of the other q.uantitles in eq.uatione (43), 

Vj, generally coinciding "with containing an arbitrary parameter. 

These simple functions reduce to the coefficients of the “natural" 
series -when, for vanish ing becomee uniiy. The arbitrary 

parameter is then so determined that the contour closes up. However, 
because of the necessity of using, to express this condition, different 
series connected "by intricate relations, this method seems to be very 

coD^licated.^ In the method studied by the author, series of the kind 
given by equations (3l)j ( 32 ), and (33) (that can be put in eln 5 )le 
relation with aeries in the exponential set) could be used, especially 
for the condition in transonic and supersonic regions where the series 
in the exponential set cease to be useful. 

It is worth'vdille to mention here that the demonstration of the 
convergence of the series in the exponential set (that Chaplygin first 
deduced in the hypergeometiic case in a somewhat complicated way) can be 
obtained very slnply and under very wide assumptions for m(w) by using 
the properties of the functions defined "by eq.uatipns (43) . Taking, for 

instance, the aeries in X^e"^^ = Immediately seen 

that Xx eauations (4-T) and (5l)> 

= (d log m/d log w)^/^ = (l — < 1 

^1 


[i - (R^l)Jd log w> 0 


and 


log 


y ow 
0 


-'-4?he problem of the closed contour has been solved In a very 
elegant way by Li^thlll (reference 8), who discovered a very simple 
development converging in all the field (subsonic and transonic) and 
gave the conditions for the closing-up of the flow behind the body. 
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Moreover j for plaiisltle m(w). 



and. 


= (1 -M^)V2 ^ 


dRj^ m 

are decreasing functions of Vj so tliat < 0 and < 0. Talcing 

now a value of n greater than 1 it is immediately seen from 
equations ( 44 ) that if_, for some value of w, E^ ^ Ej^, 

then -<^(Sn - ? , l) . ^ q gnd if 14, < E„ then >0. Hence 

if one of the two conditions is verified for some value of Wj then for 
decreasing w the value of E^ will diverge metre and more from the 

value of Ej and of E^, so that it cannot he eq.ual to unity 
at w = 0. Hence j if E^^ = 1 at w = 0, for other values of w^ 

E^ < E^_ < Ej_ and from equations (4^) 

The following limitations for are found hy use of 

eq^uatlons (47) and (52 ) , for n > l: 


e°^ <Xn< 


Naturally, these limitations hold only for real valties of E^^ and Z, 
hence for M ^1. In the same way it is proved that, for n >1, 


m°-< (Djj< e 


nZ 


The analogous demonstration for ^ and i]r^ reqxilres the assumption 
2 d(l/-w^) 

ttLat T = must te an Increasing fmetion of v. This does 

'“ 0 ® / f 0\ 

d(l/m^) 
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not seem to te too restrictive a condition, for T _^ must, in any case 

te unity at w = 0 and inf in ity at M = 1. According to tMs 
aseung)tion it can te directly verified that 


d 

mw 



_ ^ 

d(l/m^) 

dw 

, , 1/2 

1 + ^ “ ) 

\^d log w/ 

— jjL „ — . 

1 + 

\d log w/ 

2 dw 

_d(l/w^)_ 


>0 


and 


m — 


/ mw^ ^ 

m3 

^ d(l/m2) 

\m + vj 

(m. + w)^ 

d(l/w^)_ 


> 0 


By integrating ttese inequalities Between 0 and w it follows easily 
that 


m 

w 



dw^ > 



V-, 


dm.^ > 



the sign of equality holding only at w = 0. Hence equations (49) 
and (51) stow that, excluding w = 0, 1 < < T and the equation 

corresponding to (44) stows ttat like ia an increasing 

function of w. Now for -n < -1 from, tte equation in ty a 

reasoning identical to ttat developed for E-n > it can te proved 
ttat 1 < and ttat ^ > ^-1 > ^-n > 

for 1^® following Bounds told, for n > 1: 


m 




2 

m 


0 


Wj_ dnij_ 
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Similarly,, for tJtji 


X 7 ^ ^ 

The convergence of the Chaplygin (or other) series can he 
immediately demonstrated hy means of these hoods. 

For negative Integral (eventually half— integral) values of n, for 
the reasons discussed in the section entitled "The exponential set," the 
condition that E and T are uniiy at w = 0 is not sufficient to 
determine the solution. However, it is possible to give supplementary 
conditions, which are emitted for brevity, such that the resulting 
solutions of equations (^H) may he used to construct series converging 
in 1 the subsonic hodograph field exterior to a given circle. 

A new set .— Here, only briefly mentioned, is a different set of 
solutions of eq.uations (23) in closed form. It has been observed in 
the section entitled "The exponential set" that the solutions 0^ 

r\r\c\ given by equations (5O) cannot be derived from the solutions 

O 

of the eq.uations in and for n = 1. 


Conversely/ if for given G- = 0^e and G = of the 

exponential set values of F and F are deduced by means of 
eq.uation (25), or the analogous equation for the solutions obtained 
are still of the corresponding exponential set for all values of n 
but 1. 


In this case, writing instead of G = 0^® “ ^1® '*^® 

more general formulas 


G = 



WpID]- 


W3_ dm3_ + 2 


and 


G = S 


-la ( pv 


\Jvx 


“■1 ^^1 




w 


2 
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correspanding to eamtions (48), it follows from. eq.mtian, ( 25 ), 
Witt 0p = 0, that 



4m2 



i0 

+ 



are new solxrfcions of the coiTespanding eq.mtions in- ( 23 ), which are not 
included in any of the sets already discussed. 

Then if the Tariahles are inverted and interchanged and the sign 
of e is changed, it is seen that 



satisfy the corresponding eq.uations in ( 23 ); hence the respective real 
and, imaginaiy parts satisfy the equations far 0 and ■'1'’ . 

Uow applying again eq.uation ( 25 ) to the last expressions, other 
solutions of the equations in F and *¥ are found, from which hy 
changing again the variables new solutians G and. *& are found. Hence 
it is seen that the repeated application of the described process 
generates a new set of solutians F and F, and G and G. 
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Tile firsi; "tenns of "fclie corresponding Incongireasi'ble Bet of (J^j 

■wMcli can te easily obtained -by repeated application of the formula to 
vnich. eq.uatlon (25) reduces for m = w 



(■vdiere T = ve 


and successlTe Inversion of Y, are 


Y/2 



TEE GENERALIZED POTENTIAL FDHCTIONS 


In the preceding sections, it iias "been shoTm tiiat the symmetrical 
form, ottalned ty m ating the velocity and mass velocity ^and not tte 
c onne ction tetween them) appear eapllcitly in tlie Lodograph, equations, 
gives rise to an interesting general treatment of tliese eq,uations. It 
lian teen seen that the ccm^jlex functians E and G-, of which, the real 
pants are x 8^ connected hy symmetrical relations to the 

functions f* and ff, having <n and ijr as Imaginary parts. In this 
section, it will he shown that al 1 these complex functions can he 
deduced hy simple differentiations from a uniq_ue function called 
the generalized potential function. 
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A second generalized pcrtential function $ is also introdnced, 
■with, interesting properties. 


Let tlie con 5 )lex Telociiy and mass velocity te defined ty 

V = "we"^® 


> 


W = m(w)e“^ = 



j 


(55) 


and otserve that th.^ are 'bound "by the condition that their ratio 
W P 

— = — must he real and eq.ual to a prescrihed function of v, or 
V Pq 

that w must he a prescrihed function of the relative 

density* "When these conditions are satisfied, the two moduli of 
eq.uations (55) will he connected 'by IWI = m(lYl). Hence eq.tiatians (55) 
can he 'written 


so -that 


and 


=HI) 


m 
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It is now seen tliat if ¥ and Y are not 'bound "by the foresaid 
conditions and are independent, these relations can "be uned to define 
generalized complez Talues of 0(Y,¥), w(Y,¥), and m(V,¥) anfl there- 
fore of al l the related q.uantities. Hence the eq.uation 


1 _ ^ = 1 - ^ ^ = 1 - ^ ^°S (^A) 

d log w d log V d log n (¥/Y) 

defines a caii 5 )lex Mach number M = M(¥/y) idiich, if the pl^rsical 
conditions concerning Y and ¥ are satisfied, reduces to the 
real M(p/pq) and can therefore "be immediately deduced without the 
help of equation ( 56 ) "by singly replacing in the espressidn. M(p/pq) 
the real Tariahle p/Pq "by the ^nerally camplex -variahle ¥/V. 

Now the holograph eq.uatlons, caosidered for conplez values of the 
variahles, can "be tiansf armed hy taking Y and ¥ as new independent 
variahles. From eq.uations (55)j 


hv 


and 


5 m ^ dm 


1 


59 


we 


-le 


S¥ _ 0— i0 ilffi. 

5w dw 

— = e“i® 
5m 


. 5¥ _ 
S9 ■ 


-i9 


me 
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bo that 


m 


Sw \Sv dw cSW/ 


V 


^ - Y (— ^ 

Vdm. Bv 


a a B 

i — = T— +w — 

ae av aw 


y 


Hence eq.uations (23a) aad (23t) 'becanie 


(57) 




and 




•where subscripts denote partial derlTatlTes. These are satisfied If 




(58a) 




( 58 t) 


or 




^ d log w 


W 


(58c) 
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The first of these eq.xiatloriB can he satisfied if a function $ is 
introduced such that 


F 


F = *, 


and the second will also he satisfied if $ is a solution of the 
equation 

(60) 

where i-m 2 is the function of W/V defined hy eq.uation (56). 

The meaninp: of oq.uations ( 59 ) and (60) is the following. If 
a satisfying equation (60) is inown, and after and are 

CEilculated, the ri^t Talues of T and W are introduced (that is, such 
that W/V is real and IWI = m(|vl))/ then 


X = E.P.($y) 


(a = I.P.($y) 


will he solutions of eq^uations (8) . 

Observe that if in equation (60) the two ccmqjlez variables are 
replaced hy the real variables 9 and w (and m(w)) hy the inverse 
relations of equations (5.?)^ resulting equation (in a coiiq)lez ® ) 
remains unchanged hy exdhanglng w and m. The functions F and F 
(and X (d) can he deduced by relations containing the derivatives 

of 4 with respect to w, m, and 10. 
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Eq.iiations ( 23 c) and (23d) can lie treated in the same way 
result is found that if j(l/ 7 _,l/W) is a solution of 


1 



(1 - m2) 


W 


2 ^ 

2 


(Hf 


then the functions 



^ i. 

W 


will he, upon substitution of the right values of W and V, 
of eiiuations (23c) and (23d), so that the functions 




and 


t = 



will he solutions of eq.'uations (6). 


and the 

(61) 


(62) 


solutions 
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NatiiraUy 0 and ^ can also be related to From equations {.2k), 

•witli the help of equatians (57) and (58), it follows that 




& = 5wvT IT 


Henc 0 _, con^jarison with equations (62) shows that to within an 
unessential constant^ 


= V 

\V V/ 


The physical coordinates are deduced hy the following relations 
which may he obtained from equations (30) by lise of equations (55) j (57) j 
and ( 58 )* 


S + = +$ ) 

\W W YWJ 


S + IH = a-19 + 2 a J 
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After aBBlgnlng to V and Y ttelr correct yalues^ separating the real 
from, the Imaginary partj and using egiiatlons ( 2 ), It can he shoun 
ezpllcltly that 


' ° » (™w + ^ (™WW - »w«) 


which can adso he written 



By differentiating this equation, there can he deduced an 
espresslon of dz which must coincide with the expression (2), that Is 
with 


dz 




( 67 ) 


when d0 and dT|r are obtained differentiating the expression 
derived from equations (63) 



It can he shown that the eigreement exists If $ Is a solution of 
equation (60), and if W/V Is real. It seems therefore that the 
conparlson of the values of z calculated directly from equation (66), 
or deduced hy Integration of equation (67) (using equation (68)), may 
constitute an Interesting check on the accuracy of approcximate methods. 
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Cbsearve ■fclia't $ = cjY + + c^W 1 b a particular solution of 

eq.uation (6o), tlie simple meaning of -wiilcli is that constitutes an 

additive constant of od and of -4^, C 2 the same for X «-nf^ -0^ 

and c^ (cQng)lex) represents a general displacement of the origin of 
the physical coordinates. 

Particular sets of solutions of equatiotus (60) and (6l) are easily 
deduced from, the sets studied in the preceding sections hy Integrating 
the relations (59) aiid. ( 62 ) and introducing the espressions of w, 
and 9 as functions of Y and ¥. 

In the Incanpressihle case, Y = W, 1— M^=l and equations (60) 
and (6l) reduce to identities satisfied hy every function of Y. The 
■written relations reduce then to 


Fi(Y) = (Y) 


G (Y) = -Y2$i«(Y) =Y<&i**(Y) -®.'(Y) = Y^ ^ 

dY \ Y , 

zi(Y) = Pi’(Y) = %’*(Y) = 6^®(S + m) = w 


If "bhe solution corresponding to the incampressihle flow around a 
given body is kno'wn in the physical plane, ®j^(Y) can he deduced from 

these relations. The profile of the hody can, for instance, he defined 
hy a relation S = P(]l) he-tween the suhncrmal and •the normal for 1^=0 
(and hy giving this -ralue to the corresp onding streamline). 

Hence for I.P.(Ci) = 0, that is,, for I.P. (Y^j^") = I.P.($i') it 
follows that 


wS = R.P.(Y%") = wp[^I.P.(«ti') 
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Thp corresponding relation for the coii 5 )ressible case, which 
can be easily deduced from, the preceding formulas, is 

X~ ! i— 1 


for I 


p t 

M\V Wj 


= 0 , 


B-P- -I- Tv) - M 


after the correct yalues for T and ¥ are introduced. The function I 
will be the same In both cases if the profile is unchanged. 

The eq.uatl'ona (60) enfl (6l) are of the same general kind. It is 
possible to pass from, the one to the other not caHy by a substitution 
like equation (614-) but also by sin 5 >ly putting 


$ = Y¥$* 

\ ( 69 ) 


f = 


w 


f* 




since, as can be immediately verified, must satisfy equation (6l) 
nnfl eq.mtion (6o). The equations may be transformed in many ways 

by changing the two Independent variables. One of these transfcrmatlons 
is obtained by taking as new independent variables (K being a constant): 


= ^ ^ 

^ Y VJK 


(70) 
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The egiiaticm (6l) in $ (get $*) is tlms transformed into 


Tfliere the first factor of the ri^t-hand side is a function of W/V, 
hence of since 



W ^ 1 1+ (t)/|) 

T Vf 1 - (V^) 


A similar transformation may he performed on eq.uation (60) hy 
putting I = V + W an^ = T — W. 

Other interesting transformations are obtained hy taking 
5 = log Y + £(v/v) and n = 8(W/\r) and choosing in different Trays 
the functions 6 and 8. In this case the Tai*iahility of the 
T\ variahle can he restricted to the real field. Particular cases are 
obtained by taking, for instance, € = 0, 6 = log 


8 = 1 — and 8 = w(V/V) 


d log V 


( since 
\ d log 

I = X — i0 and 




An important 


(T/W) 

case is the one for Tdiich I = X — i0 and tj = X, Tdiere X is the 
same as in equation (n) or equation ( 52 ); this is obtain ed by 

. , . de ~1 dS \A - ^ 

taking , - and = r • If 


d log (J/V) (1 + Vi - M^) d log (T/W) m2 

this transformation is performed, an equation in and $ (or one 

in ® snd is deduced that may be used to obtain directly 

solutions of the kind obtained ty Bergman (reference 6) and 
Lighthill (reference 8). 

Finally, let equation (60) be TJrltten Td.th the actual isentropic 
law. In this case 

1 



+ 2L::iiM2 


a-7 
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HencOj conslcLeration of tlie obsenration following equation ( 56 ) and 
Butstitutlon of the resulting value of M(W/V) in equation (60), 
yields: 


W 








There is a corresponding equation fctr equation ( 61 ) . For real V/V" the 

p 

corresponding values of 1 — are shoTOi in figure 5 for some values 
¥ p 

of 7 "With — = — as abscissa. It is interesting to observe that 
T PQ 

if 7 is in the actual range for gasesj its value does not seem to 
affect to a great exbent the shape of the curvesj especially in the 

subsonic range. The factor 1 — is linear in ¥/V for 7 = 0 and 
in T/¥ for 7 = 2; for other values of 7 in the actual range it is 
not far from a straight line in the subsonic range. 


p 

The curve 7 = —1^ that is 1 — 11*=^ = — is also represented in 

/ 

the f igure_, and corresponds to Chaplygin* s approximation. The Kar m an — 
Tsien approximation corresponds to 1 — where the constant K 


(see equations (I 7 ) and ( 21 )) is so chosen that at infinity (P = P„) 
1 - will take the value 1 — Moo^ given by the true law. 


THE CHAITJGllJ-^MM^-H?^^ CASE 


For the Chnplygin-EAriii^--TBien approximation the ri^t-hand side 
of equation (7I) (or equation (69)) becomes zero so that the general 
solutions of equations (60) and ^,61) are 



(72) 


f = ©1(1) + © 2 (h) 
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•with. ej__, 02 j and f2 arhitraiy functicaas of the Tariahles | 
and Ti defined hy eq^uations (70) • The yalne K = 1 carrespcmds to 

Po2 

Chaply gink’s approximation; Z = (l — M , to the Karman-Tsien 


approximation.^ As it has heen ohserved in the section entitled 
"Approximate Methods" it can he convenient to chose E-values between 
the •two. 


The eoltrfcions (72) can he in^terpreted in two different •ways. 

First they can he regarded as •the exact solutions of -the corresponding 
eq.uations for a gas satisfying the ideal law (18) . 

With -the use of equations (7O) fctr | and ij, this law can he 
•written Ulhl = -H; or for real V/'T, hence for real I /tj. 


IItJ =Ih = t|T=-H (73) 

This relation allows •the expression of (and as the s^um. of two 

arhi^trary functions of | and f or of tj and 

In the second inteipretation^ equations (72) are considered as 
approximate solutions of the eq.uationB for the ac'tual law of gases. In 
this case 


U ^1 = In = = 




(7^) 


is no longer a constant^ hut is a function of w. 


^^As observed in foo^tnote 4 the actual presentation of this method 
is more coherent •than the -usual one, as the constants H and K of 
formula (18) are here deduced for a single reference condition; namely, 
the infinite point. 
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naturally tlie last espressloii, and tlie relation (T 3 )j not 
restrict tlie independenqy of | and tj; as ttey must te used only 
after the formal deductions from equations (72) of the following 
kind have teen performed.' Erom. eq.uations (72) and the application of 
eq^uations ( 59 )j 


E 


J - = (^1-f V 


^2 “ y ^2') 


and 


— + -4 = 2 fy - |fn' + 2 f 2 - Tlf 2 ’ 
T ¥ 


According to equation ( 64 ) j the last expcression coincides with 
hence, ty equations (72) j 



e2 — ^^2 ^^2 


From equations (62) , 


G = + i? = Ifi” - fi' + Tif2" - ^ 2 * 

and 

S' = *? + iT^r = ^ (^1” - - 4^2" + ^2') 
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so that 


+ im = ifi" - + Tif2" - fg' = gi(5) + ggC^i) 


and S 2 "being two functions related to ^2 


gl‘ = Sfi- 


g2‘ = ^^2’*’ 


y 


NoWj fhom eq.uation (66) for real W/V" (hence if V W 

the saBie argument -9j and | and tj the same argument 0) ; 


z = e^R 


+ le^% 


.P.< - 5fi* + I 1(1 + Tj)fi" + f 2 - 11^2' + I ^(5 ■^ 

~ ^^1* 2 ~ 'n)^’i" + ^2 “ ''1^2' ” ^ ■ 


- f 1 ~ ifjL* 2 ^^1” ■*■ ^2 2 '^^2 '*' 2 


= fl - |fi’+ \ 5%i" + f2 - Ilf2* + I ^^ 2 " + i (^ 1 ” + ^ 2 ") 


(75) 

relations 


(76) 


have 



") 


(77) 
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TMb relation gives tte physical coordinates as a function of the 

hodograph. coordinates v end 6 -idien I = I — + — 7 — 7 — 7 = 1©^® 

U m(v) >/kJ 

and Ti = ^ TT~7= 

m(v) \Jk 

law m(w) ) . 

Differentiating eq.uation (77) yields 

2 dz = + Uti|(f 2 ^"’dl + f 2 *”dTj) + (fi” + f^) dUijl 


.10 


are Introduced into it (with the assigned 


But dz is also given hy equation ( 2 ), which, may he transfctrmed in the 
following way 


2 dz = I dCJ + 1 -%= 
V ^ \f\fK 


Vk dt = I d(^^ + 1 >/f4r) + Tj d(0 - 1 s/&) 


or (see eq^uatlon (75)) 


2 dz = 5 %i'”d| + 


Hence, in accordance with the ohseivaticn following equation (68), these 
two espressions for dz coincide only if |gu| is constant; that is, 
if the law connecting I and tj is the law (73) (or (18)), hence the 
same law as the one for "vdilch the equations (60) and (61) admit the 
solutions (72) . This is what is done in the Earmanr-Tsien method. If, 
on the contraiy, the alternative ihteipretation of equation (72) is 
adopted, the coincidence ceases to ezist. In this case, if the law 
connecting 6 with n is the exact gas law, the error term, between the 
two values of z, that is, 

^ (fi" + fa") illnl 

with llql given hy equation (7^)j he regarded £is a measure of the 
approximation of idle approximate solution (72) . 
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Observe tliat if tlie law (73) Is taieiij hence if tj =-|> 

equation (75) shows that the general solntioca for 0 + i>®lr is a 
function of I alone (containing the constant H)^ that iSj it is a 
function of X — 10 ^ -where X is given hy equation (I9), a result in 
accordance -with the Chaplygin’s monogenel-ty conditions. 

The espression (77) for z can he also written hy using the 
functions and ipstead of f^ and f^ so that (see 

equation (76)), 


2z 


y + |Tll|(fl" + f2“) 




TJ dg2 + |Tilj| 






(78) 


If the law (73) is assumed, this equation reduces to 


l(igl + 




I dg - sy’ i dg 


(79) 


where g(^) = + g 2 (“S/l) is the conqjlex potential 0 + I'lSf in 

this case. In the incomqjressihle case H = 0, K = 1, nefl 2z^ =/ ^ dg 

To 0 


2 H 

with 1- = Y. Hence Tsien’s formula z = z^ — — 

immediately obtained. 



ig 

dzj 


dzj^ is 
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It is ■vrell-Jmo'wii tliat in tlie case of tlie flow round a tocly Tsien's 
focrnrula generates closed profiles only if circulation is atsent. Many 
authors have studied extentions of the method to the case with 
circulation. Bers, Germainj and leray (references 13 nni^ l4) haye 
followed a first way; Lin, Germain, and Gelhart (references 15, l4, 
and 7) a second way; here a third way of constructing flows around 
closed profiles with circulation will he shown, based on the subdivision 
of g(|) into gi(|) and g 2 (|). Let gi(TjV(zi)) = giU/Seo) "b© the 

con 5 )lex potential of an incanpressible flow around a closed profile with 
a circulation F. Then, since for | i 

(5/5J^ = 1 + k[(g/U -3 

it follows that 

dgi = ajfdzi = 0 




3^51 



U^i =r 


where the integrations are performsd along any contour in the physical 
plane enclosing the profile, or aroitnd the corresponding contour in 
the l-^lane enclosing (for sinplicity, suppose T ' and E real) . 
ITow generally, ““ 

/ (/({ - = -| r (81) 


SO that 



(82) 


and the value of z given by equation ( 79 ) for a ccmpaessible flow 
with g = gj^ is not one-valued. But if g = gj + gg far the 

given a g^ can be determined in such a way that the corresponding 

residual terms in equation ( 79 ) will compensate the value of the 
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laBt ■written integral. If -the yalue of -the cir dilation must remain 
unchanged^ must "be one-valued. TMs condition is ottained very 

sinqply ty taMng^ for instancej 


82 


= iTjy ^ 


The constant h can now he detertnlned so as to obtain •fche said 
con^jensation. From the identi-ty 


and from equation (81) it follo'ws -that 




n + r 


00 


- 1 /, 


^ = rhl„^r 


00 


Hence for real Ion and H = — I ti 

m *1 


oof 


dgg - H J = ^(^00 - Or 


If eq.iiation ( 79 ) must he one-valued, the last q.uantiiy must he 
eq[ual to (see eq.uation (82)) 
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Hence the eq.uatlons 


2tJco 

^00 


and 



will satisfy the foresaid condition of carapensaticn. 

Ohserve now that the constant factor does not depend on n. Hence 
the condition is also satisfied if (§ce/5)^ is replaced tfj P( loo/l)/P(l), 

where P is a polymonial in la/l or an infinite series converging in 
all the domain of variation of S* Hence the function 


where P is an arhitraiy function of loo/I , is analytic in all its 
domain of variation, ani^ will generate a flow around, a closed profile 
with circulation F . 

Ohseive that hy using eq.uatlons ( 70 ) for real 
expression for the constant h is 

li = -^Usa £22 vTE: _ 1 


and the 
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Jor tlie Chaplygin or E&manr-(Pslen -values of K there Is obtained^ 
respectiTely, the very simple ezpressions, both Tanishing in the 
incompressible case: 


h = 



1 


and 


h = sjl - 1 

'' 00 


Particularly simple forms of eq^uation (84) are obtained by 
putting n = 0 or n = l in eqiaatlon ( 83 ). For n = 1 an expression 
is obtained which coincides with -sdiat becomes the Lighthill solution 
for 7 = — 1 (reference 8 ) . 

The solution (84) satisfies the condition of generating solutions 
around closed profiles -v/ith circulation when jq || = Constant, and 
z is given by equation ( 79 )* H^ow if |ij || is variable, as given by 
equation (74), and the equa-bions ( 72 ) are considered as approximate 
solutions of the exact equations ( 60 ) and ( 6 I), it is sti l l possible to 
find solutions for -which z, given by equation ( 78 ), is one— valued. If 
eigaln g^( 5 A„,) conplex potential for the incompressible case, 

satisfying equations ( 80 ), -with in pleice of f, and if e^(Tf/ 7 ^„) 

is a function for which the expressions 
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(analogous to eq_uationB (80)) hold where the integrations are performed 
in the tj plane along the contour corresponding to the one of 
eq.uations (8o)j then eq^uation (82) and the ana l ogous relation for gg 

give 





2 ^ 

^00 


so that the sum 



will he zero for 



(86) 


Hence equation (78) will he one-=mlued if g 2 satisfies 
eq.uation (85) with f 2 glren hy eq.uation (86), and 


0 + i>/Et= gi(l) + g2(Tl) 
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vlll represerrtj vlien the esjiressions (7O) for I a-nfl are replaced 
"by the exact cormectloh hetween m and w, an approximate solnblon of 
the exact eq.uatlons with circulation 


r 


= r-, + r. 





z heing given hy eq.Tiatlon (78) j lli®' error term. 



representing a measiire of the approximation obtained. 

Clearly the sln5)lest way of satisfying eq^uatlons (85) and (86) is 
that of taklag^3 


82(5f) 8i(i) 


Naturally this solution still holds when hi = Constant, 

An approximate solution for the transonic case is now noted, 
corresponding to the subsonic Karmlmr-!rsien approximation. If = 1, 

then the Karm»Cn-gsien value (eq[uation (21)) for S. is zero. The 


^3if -the Karm&Or-Tslen value (eg^uatlon (2l)) is adopted for K, 
which makes the rigjht-hand side of eq^uation (71) zero at infinity, then 



00 


1 - v/l -M„2 
■ - and 

1 n/T-mJ 
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correspcmdlng cxirTe of 1 — In figure 5 redxtces to tlie horizontal 

axis 1 — 5= 0 . In. this case the ri^t— hand sides of eq^ua.tlons ( 6 o) 

artfl (6l) are zero^ anri the respective solutions are 

f = fj^(W) + Yf^Of) 

and 

5 = i ei(¥) + e^(W) 


vhere the arhitrary functions involved are hound hy the relations 


-®1 



and 


^2 




Equations (66) and (68) then give 


z = |R.P.(¥f2’) + + I.P.(f25 



^ + iTjr = -R.P.(f]^’) 


+ il.P.^ 


„ n 

V ^1 


+ ¥^2” 


+ ¥fo' - f. 


■) 


The approximate gas law corresponding to 1 — = 0 is, hy 

equation (56), m = Constant, therefore Pq/p is proportio n al to w. 

Actually m has a Tna-riminn at M = 1 , and the distance between the 
streamlines loaches there a mlnlnium. If the approximate law m = Constant 
were eidopted, the distance between the stireamlines would he unchanged 
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throu^out the field of motion, and this ■would give rise to difficulties. 
But if the solution is considered as an apiproximate solution of the 
exact eq.uations in the transonic field (even for 1 hut near 1), 

then the exact law for ia(w) can he Introduced in the solutions. This 
approximate transonic theory seems -worthy of deTelopment. 
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Figure 2.- Comparison of the . Karman-Tsien and Chaplygin approximations 

with the true isentropic law. 
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Variables a and 3 as fvtnctions of X for isentroplc gas. 
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Figure 5.- Ccaaparlson of the Chaplygin and 
the Keijmf^n-Tslen approxlmatlonfl of 
1 - as a function of p/pq vith the 
exact relation for several values of 7 . 


Figure 4. - Varlahlea and as 

functions of X]_ for 7 = 1,4, 
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